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, In this talk the results of the analytical and numerical analysis of the 

■ nonlinear Balitsky-Kovchegov equation are presented. The characteristic 

BFKL diffusion into infrared regime is suppressed by the generation of 
the saturation scale Q 3 - We identify the scaling and linear regimes for 
the solution. We also study the impact of subleading corrections onto the 
O l' nonlinear evolution. 
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1 Introduction 

One of the major challenges in QCD is the description of high energy scattering 
phenomena. In the high center-of-mass energy y/s and in the perturbative 
domain when a s <C 1 the scattering amplitude is obtained by the summation of 
diagrams leading in logs At this level, when a s is frozen, the dependence 
of the resulting cross section on the energy is governed by the power law x~ uw 
where x is the Bjorken variable. The critical exponent wp = 4 ln2a s (a s = 
a s N c /iT) is provided by the minimum of the eigenvalue function xil) °f the 
BFKL evolution kernel. 

The conceptual problem in this approach is the fact that at sufficiently high 
center-of-mass energies the BFKL amplitude violates the Froissart unitarity 
bound. This means that the validity of this approach is strongly limited and 
it has to be modified at very small x in order to guarantee the unitarity of the 
resulting cross section. 
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The solution to this problem can be provided by including the recombination 
effects, which are likely to occur at very small values of x 0. By decreasing 
the value of x at fixed gluon virtuality k%, the density of partons becomes so 
large that they start to overlap. In this case the gluon splitting process must be 
supplemented by a competing gluon recombination. In terms of evolution the 
master equations become nonlinear with an additional quadratic term which 
suppress the growth of the amplitude with energy and restore the unitarityQ. 

There have been extensive studies on this problem, see J|, @-|14 , which 
result in the nonlinear evolution. One of the important outcomes of these studies 
is the existence of the saturation scale Q s (x) which is a characteristic scale at 
which the parton recombination effects become important. In particular case of 
the Balitsky-Kovchegov equation [|l4| the existence of such scale has yet another 
important impact on the picture of the BFKL evolution. The diffusion into the 
infra-red, which is the charactersitic property of BFKL evolution is strongly 
limited due to the existence of the saturation scale |j| |23|. In fact, in the regime 
when the gluon transverse momenta k < Q s (x), the solution to the nonlinear 
equation Q becomes a function of only one combined variable k/Q s (x) 
[ p4[ . In the regime of high momenta, k > Q s (x) the parton density is small and 
the evolution is governed by a normal linear equation. 

In this talk we present the analytical and numerical analysis of the Balitsky 
- Kovchegov |l4| equation which is a nonlinear evolution equation in the leading 
logs limit. We illustrate the emergence of the saturation scale and scaling and 
show that it leads to the suppression of the infra-red diffusion. We also consider 
the case with additional NLL effects such as kinematical constraint and running 
coupling. 

The results presented in this talk have been obtained in the collaboration 
with K. Golec - Biernat and L. Motyka. For the details of the calculation the 
reader is referred to 031. 



2 Nonlinear evolution equation 

The Balitsky-Kovchegov equation (l4) has been derived as an evolution equa- 
tion for the dipole-nucleus amplitude in the dipole picture by a summation of 
multiple Pomeron exchanges in the leading logs level and in the large N c limit. 
The resulting evolution equation reads 



aV(r,b,Y) = _ (Jfg)JV)(rib>y) _ 



dY 

•72 



/rj^v 1 r v l v A- v 1 

^^r^i^ + ^b + ^n^b + ^r), (i) 

2 Recently |3| it has been pointed out that th e si tuation can be actually more complicated 
in a sense that the Balitsky-Kovchegov equation |14[ could lead to the local saturation but not 
to the unitarisation due to the fact that the target radius in impact parameter space could 
grow as fast as a power with energy. 
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where a s = N c a s /-K, and the linear term is determined by the BFKL kernel 

(K®N)(r,b,Y) = 

f d 2 r' ( r 2 r' r 2 1 

= / -1^ 7 T^N(r + r',h+-,Y)- — 9 - -=N{r,b,Y)}. (2) 

J nr' 2 \(r + r') 2 y 2 ' r' 2 + (r + r') 2 J 

The function N(r, b, Y) is the dipole-nucleus amplitude for the scattering of the 
dipole with transverse size r at impact parameter b and at rapidity Y. 

In the linear approximation, when each dipole scatters only once off the 
nucleus, the BFKL equation in the dipole picture is obtained. The non-linear 
term in (Q) takes into account multiple scatterings and is essentially determined 
by the triple pomeron vertex |l3[ in the large N c limit. Eq. (^) unitarizes the 
BFKL pomeron in the sense that at x — > and Q 2 fixed, 

F 2 ~ Q 2 hx(l/x) . (3) 

Thus the power-like rise with energy for the BFKL pomeron is tamed jl5| . 

For the subsequent analysis we shall assume the approximation of the big 
nucleus, i.e. r <C b which allows us to factorize the impact parameter depen- 
dence in Eq. (jl|) . We also consider the spherical symmetric solutions in r and 
transform the equation ([j]) into the momentum space by performing the Fourier 
transform 



; d 2 r N(r Y) f°° dr 

cj>{k,Y) = I ^-exp(-ik-r)^-i = / - J Q (kr) N(r,Y), (4) 
z7r r ,/n r 



where Jo is the Bessel function. In this case the following equation is obtained 

= a s {K>®<t>){k>Y) ~ a^ 2 (fc,F), (5) 

and the action of the BFKL kernel is given by 

dk 12 f k 12 4>{k' \Y) - k 2 4>{k,Y) k 2 4>{k,Y) 



poo jy2 ( 



\k 2 - k' 2 \ V4fc' 4 + fc 4 



(6) 



where now k and k' are the transverse momenta of the exchanged gluons in the 
BFKL ladder. 



3 Saturation scale and geometric scaling 

In order to study the diffusion and scaling properties of the solution we shall 
consider the function k<j>(k,Y). In the case of the linear BFKL equation this 
function is a Gaussian concentrated around some initial scale fco and with width 
increasing with rapidity, leading to a diffusion. In Fig.[l] we illustrate this dis- 
tribution for the case of the solution to the linear BFKL and in the nonlinear 
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Figure 1: The functions k<p(k,Y) constructed from solutions to the BFKL and 
the Balitsky-Kovchegov equations with the delta-like input for different values 
of the evolution parameter Y — ln(l/x) ranging from 1 to 10. The coupling 
constant a s = 0.2. 



Balitsky-Kovchegov equation for different choices of rapidity. As an initial con- 
dition we have chosen a delta function 5(k— ko). In the case of the linear BFKL 
evolution the solution is always peaked at k = ko and exhibits the well known 
diffusion pattern, in which the momentum distribution has an increasing width 
with growing rapidity. In the nonlinear case however the solution behavies quite 
differently. The peak of the distribution fc max moves towards higher k values 
as the rapidity increases and the solution becomes washed out from the k < ko 
regime. Only in the region k ^> fc max it coincides with the linear evolution. 

The impact of unitarization of the BFKL pomeron on the infra-red behaviour 
can be also visualised by studying the properties of the following normalised 
distribution 

m y) = k ^ Y) (7) 

and by performing the projection of this function onto the (hik/ko,Y) plane, 
Fig. | 
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Figure 2: The lines of constant values of the BFKL and the BK re-normalized 
solutions ^>(k,Y) (Y = \n{l/x)) in the (log 10 (fc), \og 10 (l/x)) -plane. 



Again, for small Y , when the non-linearity in the BK equation is negligible, 
the re-normalized solutions (Q) of the BFKL and the BK equations coincide. 
With increasing Y, when the non-linear effects become important, the difference 
between them in the region of small k becomes fully visible. Note that in certain 
region of (lnfc/feo, Y) in Fig. || the contours become parallel straight lines. This 
means that ^(k,Y) in this region is a function of the combination 



£ = ln(fc/fc ) — AF = In 



k 



fco exp(AF) 



A > 



(8) 



instead of k and Y separately. This is referred as the geometric scaling property. 
We note, that this scaling holds in the regime where k < k max and becomes 
violated when k > fc max . This suggests that we can identify the saturation scale 
Q S (Y) with 

n max 



3 (Y) = k max (Y) = Q cxp(Ar) 



(9) 



with the exponential dependence on rapidity governed by the value of the scaling 
parameter A. 
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Figure 3: The function (k/Q s (Y)) <f>(k, Y) plotted versus k/Q s {Y) for differ- 
ent values of rapidity Y ranging from 10 to 23. The saturation scale Q S (Y) 
corresponds to the position of the maximum of the function k(f>(k, Y). 

The solution to Eq. (^) has the same scaling property as the function ip(k, Y) 
namely 

<t>(k,Y) = <f>(k/Q a (Y)) , (10) 

provided that <j>{k max ,Y) = const. We have checked that this condition is 
satisfied. 

We have checked that the scaling coefficient A (0) is a universal quantity and 
it does not depend on the type of the input distribution whereas the normali- 
sation Qq is determined by the initial conditions. 

The transition between the scaling and linear regime can be perhaps best 
seen in Fig. || where the function k/Q s (Y) (f>(k,Y) is plotted as a function of 
the scaling variable k/Q s (Y) for different values of the rapidity Y. The scaling 
behaviour is represented by a common line to the left of the point k/Q s (Y) = 1. 
The slow departure from scaling for k > Q S (Y) is clearly visible. In this sense 
the line k = Q S (Y) is only an approximation characterizing the position of the 
transition region in the (fc, Y~)-plane. However, the choice based on the position 
in k of the maximum of k(f>(k, Y) as a function of Y is the most natural one. 
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From the numerical solution we have extracted the value of the scaling pa- 
rameter and found A ~ 2.05 a s which is in agreement with the estimates of 
[18, [H| for the Balitsky-Kovchegov equation and also with the previous work 

The geometric scaling is also the property of the Golec-Biernat and Wusthoff 
saturation model |25| which succesfully described the data on inclusive and 
diffractive proton structure function. In this model it is assumed that the whole 
energy dependence of the dipole cross section cr^r, Y) = a d{r 2 Q 2 S {Y)) is driven 
through the saturation scale Q S {Y) (or saturation radius 1/Q S (Y)). It then 
results in an approximate scaling property of the data for the total photon- 
proton cross section 

a rp {x,Q 2 ) = a Tp (Q 2 /Q 2 s {Y)). (11) 
Such scaling law was found in the small- x DIS data p(| . 



4 Analysis beyond LL level 

The Balitsky-Kovchegov equation has been derived in the leading logs level, 
therefore it is important to study the impact of the subleading corrections. 
Although it has been argued |[7j that the unitarity corrections might become 
important before next-to-leading ones, the study of the latter ones is crucial due 
to the large numerical value of the subleading series. After the evaluation of the 
NLL contribution to the BFKL amplitude has been performe d p8[ it turned 
out that the correction is very large and need to be resummed ]29| in order to 
stablize the final result. 

We have tested two important types of the NLL corrections: the inclusion 
of running of the coupling a s (k 2 ) and the so called kinematical constraint. 

In the case of linear BFKL the running of the coupling poses serious prob- 
lems due to the existence of the Landau pole and thus it is necessary to reg- 
ularise a s (k 2 ) at small scales. This results in a strong dependence of the so- 
lution on the cut-off (or freezing ) parameter fco- The intercept of the BFKL 
Pomeron turns out to be dominated by the values of a s {k'^ > ) and instead of 
the typical diffusion pattern one has a factorized behaviour of the solution 
k(j}(k,Y) ~ exp{Ay} i [In (fc 2 /^)]" at large rapidities. The distribution of 
the gluon momenta is dominanted by the virtualities in the infrared regime 
~ ko- 
lw the case of the Balitsky-Kovchegov equation with the running coupling one 
has still problem of regularisation around the Landau pole, however the solution 
is itself much more stable with respect to the details of the phenomenological 
regularisation. It turns out, that since the saturation effects are very strong 
in the regime of small values of k, they tend to decrease the rapid rise of the 
amplitude for the values of running coupling evaluated at scales around the 
cutoff scale fco- The rapidity - dependent saturation scale, Q S (X) shifts the 
momentum distribution out of the infrared regime into the perturbative domain. 
This phenomenon can be best visualized by means of similar contour plots as 
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Figure 4: The re-normalized solution ^(fcjV) for the Balitsky-Kovchegov equa- 
tion with the running coupling constant and the infra-red cut-off k$ = 0.1 GeV 2 
as function of log w (l / x) and log 10 (fc/l GeV). 



before, Fig.^. In the linear BFKL case one observes that the solution initially 
undergoes the diffusion pattern, slightly asymmetric due to corrections coming 
from the running coupling, and then suddenly it drops into the small scales 
regime k ~ fco where it exhibits factorized (in Y and Ink /kg) behaviour. On 
the contrary, the Balitsky -Kovchegov equation which initially coincides with 
BFKL, also undergoes some form of transition but then the distribution moves 
away from the infrared regime due to the generation of the saturation scale. We 
have estimated the saturation scale in the case of Balitsky-Kovchegov equation 
with running coupling. By taking the ansatz that the local exponent X(Y) = 
dln(Q a (Y)/A) /dY takes the form X(Y) = 2a s {Q 2 s (Y)) where A = A QC d we 
have found that the saturation scale is given by the solution to the differential 
equation 

dHQs(Y)/A) _ 12 

dY & ln(Q s (V)/A) ' 1 ' 

with the initial condition Q s (Yq) = Qq and Yq chosen in the region where scaling 
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Figure 5: The function (k/Q s (Y)) <p(k,Y) in the running coupling case, plotted 
versus k/Q s (Y) for different values of rapidity Y ranging from 15 to 32. The 
saturation scale Q S (Y) is taken from equation with the initial condition 
Q S {Y = 0) = 2.0GeV. 
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sets in. The solution takes the form 



Q s (y)=Aexp( J^(Y-Y )+Ll) , Y > Y , (L'Si 



where Lq = ln(Qo/A). It follows that the local exponent X(Y) decreases with 
increasing rapidity, and \(Y) ~ 1/VY for very large Y. 

Such dependence is indeed seen in the numerical analysis. 

In Fig. H we illustrate scaling in the case with running coupling by showing 
the function (k/Q s (Y)) cj)(k,Y) plotted versus k/Q s (Y) for different values of 
rapidity. Q S (Y) is given by formula ( fl3| ) with the initial condition Q S (Y = 
0) = 2 GcV. The overlapping curves at low values of the scaling variable clearly 
indicate that for k < Q S (Y) scaling is satisfied to a very good accuracy, thus 



justifying our ansatz fll3| ) for the saturation scale. We have also tested the 
impact of the so called kinematical constraint onto the solution of Balitsky- 
Kovchegov equation and found that though qualitative features are the same 
(with generation of the saturation scale and scaling) the absolute numerical 
value of the solution is strongly decreased and also the saturation scale is shifted 
towards smaller value of gluon momenta. 

To summarize, the diffusion into infrared in the case of the Balitsky-Kovchegov 
equation is strongly suppressed due to the emergence of the saturation scale 
which can be approximately identified with the maximum of the momentum 
distribution k</)(k,Y). The solution to the nonlinear equation has the property 
of the geometric scaling in the regime where k < Q S (Y) whereas in the case 
when k > Q S (Y) the solution enters the linear regime. In the case of running 
coupling the nonlinear equation turns out to be more stable as compared with 
the linear BFKL evolution. The sensitivity to the treatment of the infra-red 
region is much smaller than in the linear case due to the appearance of the 
saturation scale. In contrast to the BFKL equation with the running coupling, 
the large Y asymptotics of the dipole scattering amplitude arising from the BK 
equation is governed by gluon momenta in the perturbative domain. 
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